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FOREWORD 


The  number  of  publications  in  the  field  of  pressure  vessel  stress 
analysis  has  increased  precipitously  in  recent  years  and  the  information 
contained  therein,  though  pertinent,  is,  many  times,  lost  in  this  sea  of 
literature.  In  addition,  a  large  proportion  of  the  solutions,  as  presented, 
are  not  carried  to  the  point  where  they  can  be  used  directly  for  determining 
stress  distributions.  The  purpose  of  this  report,  then,  is  to  make  avail¬ 
able  a  compact  and  adequate  summary  of  the  formulas  and  to  present  them 
in  a  readily  usable  fashion  for  the  stress  analysis  of  pressure  vessels  that 
are  commonly  encountered  in  missile  design.  That  such  a  report  as  this 
derives  largely  from  the  work  of  others  is  self-evident,  and  it  is  the 
author's  hope  that  due  acknowledgement  has  been  made  of  the  immediate 
sources  of  all  material  here  presented.  However,  while  most  publications 
deal  with  head  closures  joined  to  long  circular  cylindrical  shell  sections, 
this  report  places  equal  emphasis  on  designs  in  which  the  circular 
cylindrical  shell  section  is  short  and  is  integrally  joined  to  head  closures 
of  different  shapes  and  thicknesses.  The  equations  are  organized  such 
that  interchangeability  of  head  closures  is  a  matter  of  algebraic  manipu¬ 
lation.  It  is  believed  that  the  material  contained  in  this  report  will  be 
helpful  to  those  involved  in  pressure  vessel  stress  analysis  and  will  be 
particularly  useful  for  members  of  the  Structures  Section,  Solid  Mechanics 
Department,  who  are  often  asked  to  appraise  pressure  vessel  design 
speedily  and  yet  with  sufficient  accuracy  as  to  insure  a  degree  of  confidence 
in  their  appraisal.  The  material  presented  herein  can  also  be  used  to  spot 
check  stress  levels  obtained  from  computer  solutions. 

The  author  would  like  to  take  this  opportunity  to  express  his  thanks  to 
Mr.  Walter  Smotrys  of  the  Structures  Section  for  his  assistance  in  checking 
the  algebraic  manipulations  associated  with  the  preparation  of  this  report. 
He  also  wished  to  thank  Miss  M.  J.  McNeil  for  her  computations  of  the 
frequently  encountered  parameters  and  the  preparation  of  the  associated 
graphs. 
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Finally,  it  should  be  said  that,  although  every  care  has  been  taken  to 
avoid  errors,  it  would  be  over -sanguine  to  hope  that  none  had  escaped 
detection;  for  any  suggestions  that  readers  may  make  concerning  needed 
correction,  the  author  will  be  grateful. 
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STRESSES  IN  THIN  VESSELS  UNDER  INTERNAL  PRESSURE* 

ABSTRACT 

Elastic  stresses  in  thin  shells  of  revolution  under  the  action  of  internal 
pressure  are  presented.  The  formulas  given  are  developed  on  the  basis 
of  Love's  classical  shell  theory.  The  pressure  vessel  configurations 
under  consideration  consist  ofvarious  commonly  encountered  head  clo¬ 
sure  designs  integrally  joined  to  circular  cylindrical  shell  sections  which 
may  be  classified  as  long  where  the  characteristic  length  of  (3cf  >  4.  0  or 
short  where  p^f  <4.0.  In  addition  to  the  membrane  stresses,  the  bending 
stresses  resulting  from  forces  and  moments  at  the  junctures  of  the  heads 
and  cylinders  are  also  presented.  The  concept  of  edge  influence  numbers 
is  used  where  convenient  to  express  the  discontinuity  forces  and  mo¬ 
ments  at  the  junction.  Many  important  parameters  are  expressed  in 
graphical  forms  to  facilitate  analysis. 
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1.  INTRODUCTION 


In  thii  report,  &  pressure  vessel  will  be  defined  as  a  container  that 
must  withstand  an  internal  pressure,  and  is  in  the  form  of  a  surface  of 
revolution  with  a  wall  thickness  small  compared  to  the  radii  of  curvature 
of  the  wall  (h/a  £  1/15).  It  is  well  known  that  in  every  boost  vehicle  for 
space  launching  and  ballistic  missile  systems,  the  pressure  vessels 
comprise  a  large  percentage  of  the  total  system  structural  weight.  Con¬ 
sequently,  their  design  merits  considerable  attention  from  the  analyst.  In 
flight,  liquid  systems  will  be  pressurized  either  by  the  acceleration  forces 
acting  on  the  hydrostatic  head  of  the  propellant  plus  a  pressure  head  for 
efficient  pumping  of  the  propellants  to  the  engines  or,  in  the  case  of  solid 
propellants,  by  the  pressure  of  the  burning  propellant  gases.  In  this 
report,  however,  attention  will  be  confined  to  the  effects  of  loadings  caused 
by  internal  pressure  only. 

Pressure  vessels  with  very  thin  walls  which  offer  no  resistance  to 
bending  would  be  subjected  only  to  direct  stresses  uniformly  distributed 
through  the  thickness.  In  other  words,  the  wall  is  acting  as  a  membrane. 

The  associated  stresses  are  called  "membrane  stresses." 

However,  when  the  wall  offers  resistance  to  bending,  bending  stresses 
occur  in  addition  to  the  membrane  stresses.  In  pressure  vessel  design, 
bending  stresses  arise  as  a  result  of  (a)  change  in  curvature,  (b)  change  in 
slope,  and  (c)  change  in  wall  thickness.  Bending  stresses  resulting  from 
these  causes  are  called  "discontinuity  stresses."  These  stresses  are 
obviously  most  severe  at  or  near  the  discontinuity.  They  do  not  vary 
circumferentially  because  of  the  presumably  axial  symmetry  of  the 
structures  and  decay  rapidly  to  negligible  values  in  a  distance  I  >  4p 
measured  from  the  point  of  discontinuity  on  the  meridional  arc. 

From  the  above  brief  description,  it  is  clear  that  in  order  to  determine 
the  complete  state  of  stress  in  a  pressure  vessel,  it  is  necessary  to  find 
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the  membrane  and  the  discontinuity  stresses.  These  separate  effects  may 
be  linearly  superimposed. 

It  is  well  known  to  those  familiar  with  the  calculation  of  bending 
stresses  in  shells  that  the  process  involved  is  both  complex  and  time 
consuming.  To  a  great  extent,  the  complexity  is  associated  with  the 
determination  of  the  discontinuity  circumferential  shear  forces  and  bend¬ 
ing  moments  existent  at  the  juncture  of  head  closures  and  cylinder  or  at 
other  discontinuities  mentioned  earlier.  Accordingly,  a  large  portion  of 
this  report  will  be  devoted  to  the  calculation  of  these  discontinuity  forces 
and  bending  moments.  The  concept  of  edge  influence  numbers  is  used 
where  convenient  for  this  purpose.  (By  "edge  influence  number1’  is  meant 
the  displacement  or  rotation  of  the  edge  of  the  shell  due  to  unit  values  of 
the  edge  bending  loads  or  unit  values  of  the  pressure.)  This  concept  is  not 
new.  Indeed,  Muller- Breslau  used  it  many  years  ago  for  beams,  trusses, 
etc, ,  and  other  authors  have  extended  it  to  shells,  such  as  the  works  of 
Galletly  (Ref.  1);  Watts  and  Lang  (Ref.  2);  Taylor  and  Wenk  (Ref.  3);  to 
mention  a  few.  Having  the  edge  influence  numbers  will  greatly  simplify 
the  stress  analyst's  task  in  formulating  the  compatibility  equations  at  the 
junctions  of  head  closures  and  cylinders. 

Appropriate  expressions  have  been  included  in  this  report  to  permit 
the  determination  of  stress  distribution  throughout  the  shell. 

The  following  assumptions  are  employed  throughout  this  report; 

(a)  The  material  under  consideration  is  assumed  to  be  a  perfectly 
elastic,  homogeneous,  and  isotropic  solid. 

(b)  The  load  is  assumed  to  be  entirely  due  to  internal  pressure,  so 
that  support,  dead  weight,  and  similar  loads  are  completely 
neglected  in  this  report. 

(c)  All  heads  are  considered  to  be  complete  (without  holes)  and  to  be 
free  from  any  stress  raisers  other  than  the  head- cylinder 
juncture  itself. 
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(d)  Each  cylinder  and  head  closure  is  assumed  to  have  uniform 
thickness  (although  the  thicknesses  need  not  be  the  same). 

(e)  The  middle  surfaces  of  the  cylinder  and  head  at  the  juncture 
are  assumed  to  be  continuous. 

Before  we  proceed  to  the  main  portion  of  the  report,  it  is  important 
to  recognize  that  the  formulas  presented  are  derived  on  the  basis  of  the 
elastic  behavior  of  the  shell  materials,  so  that  when  the  highly  localized 
discontinuity  stresses  calculated  by  these  expressions  exceed  the  elastic 
limit  of  the  material,  local  plastic  deformation  in  the  region  of  the  junction 
of  the  cylinder  and  head  will  take  place.  This  plastic  deformation  prevents 
stresses  of  the  calculated  intensities  and,  in  areas  where  it  is  possible  for 
plastic  flow  to  occur,  the  maximum  stress  may  be  limited;  the  high  peaks 
indicated  by  elastic  analysis  will  generally  be  redistributed  and  increase 
the  stress  at  nearby  points. 

One  might  therefore  challenge  the  value  of  the  elastic  analysis.  The 
following  statements  may  be  made  in  its  defense: 

(a)  The  high  stresses  indicated  by  elastic  analysis  are  valuable 
since  they  point  out  potential  trouble  spots. 

(b)  The  elastic  analysis  is  needed  to  define  the  areas  where  plastic 
flow  may  occur. 

(c)  Although  high  local  stresses  in  steel  or  aluminum  vessels  are 
frequently  relieved  by  plastic  flow,  a  good  design  cannot  always 
rely  on  ductility  as  insurance  against  the  bad  effects  of  sharp 
corners  and  other  stress  raisers.  For  low  temperature 
operations,  some  normally  ductile  materials  become  brittle 
and  failure  as  a  result  of  this  brittleness  can  occur. 
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2.  MEMBRANE  THEORY  SOLUTIONS 


A  tabulation  of  Membrane  Theory  eolutione  for  preaeure  vessel 
configurations  of  practical  interest  under  the  action  of  internal  pressure  is 
offered  in  this  section. 

2.1  Thin  Circular  Cylindrical  Shell  Under  Internal  Pressure 


ttr‘%-(l‘l) 

where  ^  =  hoop  membrane  stress 

=  meridional  membrane  stress 
=  radial  displacement 
p  =  uniform  internal  pressure 


a  a  mean  radius  of  curvature 
h  a  shell  thickness  „ 

E  >  modulus  of  elasticity 
v  a  Pois.von's  ratio 

2.2  Thin  Spherical  Shell  Under  Internal  Pressure 


ffe  s  %  =  2h 

where  =  hoop  membrane  stress 

=  meridional  membrane  stress 
=  radial  displacement 
p  =  uniform  internal  pressure 


a  »  mean  radiue  of  the  apherical  ahell 
h  »  eheU  thicknen 
E  »  modulo*  of  elasticity 
v  ■  Poieeon'a  ratio 

2. 3  Thin  Ellipeoidal  Shell  Under  Internal  Preeeure 


where  Cg  =  hoop  membrane  etreae 

*  meridional  membrane  etreea 
uf  *  radial  dieplacement 
p  >  uniform  Internal  preeeure 

a  «  eemi-major  axie  of  the  ellipaoidal  head,  meaeured  to  the 
middle  eurface  of  the  ahell 

b  *  eemi-mtnor  axie  of  the  ellipeoldal  head,  meaeured  to  die 
middle  eurface  of  the  ehell 

h  -  ehell  thickneee 

E  s  modulua  of  elaeticity 

v  =  Poieaon'a  ratio 

*«| 

^  I  Principal  radii  of  curvature  of  the  ellipaoid 

X1 

I  Carteaian  coordinate  ayatem 

y  1 

2. 4  Thin  Conical  Shell  Under  Internal  Preaeure 
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ffe  s  *r tan  ° 

sStano 


where  <r^  *  hoop  membrane  stress 

o,i  *  meridional  membrane  stress 
uy  =  radial  displacement 
p  =  uniform  internal  pressure 
h  =  shell  thickness 

s  =  distance  of  a  point  of  the  middle  surface  from  the  vertex 
measured  along  a  generator 

E  =  modulus  of  elasticity 

v  £  Poisson's  ratio 
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K  . 

\  * 
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O  2. 5  Thin  Elastic  Simply  Supported  CtrcuUr  Plate  Undar  Uniform  Load 

Over  Entire  Surface 


? 


At  R: 


L, 


*.*£(*♦4-4) 

T  8h2  '  a  / 

3.,^] 

y  =  3P.a2l\.-.yj[4  (5  + v)f  (i  +  v)  -  R2  (J  +  v) 

8Eh3  l  2  2a2 
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At  center:  a-*,  and  y  attain  their  respective  maximum  values,  and  are 

■  — "  T  “ 

given  by: 


'r  s*9 
max  max 


_  3pa4  (1  -  v)  (5  4  v ) 

max'  16  Eh3 


where  <r  =  unit  stress  at  surface  of  plate  in  the  radial  direction  and  is 

r  positive  with  tension  at  the  lower  surface  and  equal  compression 
at  the  upper  surface 

0*0  =  unit  stress  at  surface  of  plate  in  the  tangential  direction  and  is 
positive  with  tension  at  the  lower  surface  and  equal  compression 
at  the  upper  surface 

y  a  vertical  deflection  of  plate  from  original  position  and  is  positive 
for  downward  deflection 

p.=  uniform  pressure  o\'er  the  surface  of  the  plate 
h  =  plate  thickness 
E  =  modulus  of  elasticity 
v  =  Poisson's  ratio 
a  =  outside  radius  of  plate 

R  a  distance  to  any  given  point  on  surface  of  plate 
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3.  BENDING  THEORY  SOLUTIONS 


A  tabulation  ol  B«nainl  Theory  Solution,  lor  prowur.  «•••> 

01  practical  intact  «*••'  th.  action  .1  varhc  a<|.  1«* 

is  offered  in  thie  section. 

3  1  Circular  Cylindrical  ShejU 

Th.  following  .l,«.  convention  for  th.  roUtlon  and  d.n.c«on  *'«  »“* 

Positive  Deflection  -  radially  inward  with  reepect  to  the 
Positive  Deue  center  line  of  the  cylinder 


Poiitive  Rotation 
These  signs  are  graphically  represented  below: 


clockwise  viewing  the  upper  cut  of 
cylinder 


The  following  nomenclature  is  used  in  this  section: 

*  s  meridional  stress 
x 

ff_  =  circumferential  stress 
*0 

a  =  radial  displacement 
r 

0  =  slope 
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ft4  -  W'.* ) 

p  -  rr 

a  h 

D  -  flexural  rigidity  of  the  shell 

_  Eh3 
12(l-vZ) 

Where  ±  signs  occur  in  the  expressions  for  stresses*  the  upper  sign  refers  to 
the  inner  surface  of  the  cylinder  and  the  lower  sign  to  the  outer  cylinder. 

The  following  notations  are  introduced  for  convenience: 

Notation 


J 


n 


n. 


n. 


ne 


n. 


nB 


‘10 


a 


n 


Description 

Plotted  on 

e*^x  sin  px 

63 

e~^X  cos  px 

64 

"^X(sin  px  +  cos  px) 

65 

"^X(cos  px  -  sin  px) 

66 

sin  px  sinh  px 

67 

sin  px  cosh  px 

68 

cos  px  sinh  px 

69 

cos  px  cosh  px 

70 

sin  sinh 

71 

sin  pi  4  sink  pi 

sin  cosh  QL 

sin  pi  +  sinh  pi 

72 

cos  sinh 

sin  pi  4  sinh  pi 

73 

J 


J 
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NoUtlon 

Description 

Plotted  on  Page 

cos  ^  cosh  ^ 

74 

°12 

.Tn’pr+TliKFI 

°13 

linh  pi  -  tin  pi 
•in  pi  -I-  •ink  pi 

75 

°14 

coah  Si  +  coa  pi 
•in  pi  +  ainK  pi 

76 

°15 

coah  pi  -  coa  pi 
•in  pi  +  ainK  Pi 

77 

°16 

r 

sinh  pi  cosh  pi  -  sin  pi  cos  pi 

78 

sinh^  pi  -  sin*  pi 

°17 

sinh  pi  cosh  pi  +  sin  pi  cos  pi 
ainh^  pi  -  ain^  pi 

79 

°18 

coa  pi  ainh  pi  -  sin  pi  coah  pi 

80 

1  2 
sinh  pi  -  sin  pi 

fl19 

cos  pi  sinh  pi  +  sin  Bi  cosh  pi 

81 

sinh*  pi  -  sin  pi 

n? 

o 

sinh^  pi  +  sin*  pi 
sinh6  pi  *  sin  pi 

82 

°2I 

sin  pi  sinh  pi 

- 2  "  2 — 

sinh  pi  -  sin  pi 

83 

«22 

•inZ  pi 

- . — — c - - - 

sinh6  pi  -  sin6  pi 

84 

°23 

sinh*  pi 

- -  ■■■  ■■■■  - £ — 

sinh6  pi  -  sin  pi 

85 
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3. 1.  b  Uniform  Radial  Moment  M_  In-Lb  Per  Linear  Inch  of  Circumference 
■ 11  ■■ 0  . .  —  — 1  ■■  —  "  1111  ■  111 

at  End  of  Long  Circular  Cylindrical  Shell 


Mq 

0=  pC-n2 


r(x=o) 


M. 


20fcD 


Mo 

V=o)  =  pF 


2)  Unequal  Edge  Forcee 

a, 

°u 

rr 


-  A,.a.)As] 

2&  OL 

(Jl,l,a,+Jl.,tal)A,] 

2?  +CA„a.*A,a,XA,-At)] 
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ur  (x=0)  =  ‘  ^3^  (nl8Q2  +  °16Q1) 
°(x=0)  '  ’  (2fl21Q2  *  aZOQl] 

Ur(x=<)  =  *  “Td  <°16Q2  +  «18Q1> 
®(x=l)  =  *  (Q20Q2  '  2fl21Ql) 


3. 1.  d  Short  Circular  Cylindrical  Shell  Bent  by  Moments  Distributed 
Along  the  Edges 

1 )  Equal  Edge  Moments 


Mo 


M0 
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r,  “  ^  {?  a[jl,,)(Jlt+Ai)+  -a,.(%A»)  i 

■xKC^t)-il„(A6-A,)]} 

2|i 

i  m.a,s 

2)  Unequal  Edge  Moments 

M,  Mt 


l 
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J 


^  h  |^(^i  ■^’17“  Mjll  raX‘^,ft"‘^’7)~-^’^^i'^'x*‘  + 

~Ui‘~~  *|J*0  ^ ^ ,^*1 17_  ^x‘^’  »)C'^’4+  ^  7)  ~  (^  I  ^»"  2W*  ^•il)-^,g“  ^1*^5  j 

-  m.CV-M] 

=  “  ip*^[  ^i^***”  2KAi-^*n] 

^ (>•■»)  po  C^i'^17”  ^z*®*o) 

^(2M« V  M A.) 

®(*-£)  =  “~pD  ^2-^17) 

3.  2  Shells  of  Surfaces  of  Revolution 

The  following  signs  convention  for  the  rotation  and  deflection  are  used: 

Positive  Deflection  -  radially  inward  with  respect  to  the  axis 
of  revolution  of  the  shell 

Positive  Rotation  -  clockwise  viewing  the  upper  cut  of  the 
shell 

j 
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The  nomenclature  used  in  this  eection  it: 


o\  =  meridional  iitreti 
o'q  =  circumferential  ttrets 
=  radial  displacement 
0  =  slope 

0  =  angle  between  axis  of  revolution  and  normal  to  wall 

R  =  mean  radius  of  circumference  of  discontinuity  circle 

=  radius  of  curvature  of  the  section  perpendicular  to  the  meridian 
at  the  point  in  question 

=  radius  of  curvature  of  the  meridional  section 


c4_  3(  1-v  ) 
P - TTT 


D  =  flexmal  rigidity  of  the  shell 

_  Eh3 
"  12(1- v^) 

Q  =  uniformly  distributed  circumferential  shearing  force  at  the 
°  shell  edge  perpendicular  to  the  meridian 

M  =  uniformly  distributed  circumferential  bending  moment  at  the 
shell  edge 

t  =  distance  measured  along  the  meridian  of  the  shell  from  the  edge 
to  the  point  in  question  and  is  positive  away  from  the  junction. 

£  =  summation 


v  =  Poisson*s  ratio 


Notations  in  Shell  Analysis 

Where  *  signs  occur  in  the  expressions  {or  stresses,  the  upper  sign  refers 
to  the  inner  surface  of  the  shell  and  the  lower  sign  to  the  outer  surface  of  the 
shell. 

3.  2.  a  Uniform  Normal  Shear  Q  Lb  Per  Linear  Inch  of  Circumference 

■  .  o  ■  ...  -  —  — . .  .i  .  - n  —i  -n  -  ■ 

at  End 
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3.  2.  b  Uniform  Radial  Moment  M  In -Lb  Per  Linear  Inch  of  Circumference 


6RM 

*  n  R,  3 


2M 

ff0  =  TC^IT 


ZR2n4  ±  £  vro3 

► 


(1-V  )cot  fe 
— * — 3 


4 


M 

o 

u  At  - *■—  sin  0 

r|‘s°>  2pD  ° 


M 


,(6=oj =  ’  p“ 


In  applying  the  above  formulas  for  the  determination  of  stresses  and 
displacements  for  shells  of  surfaces  of  revolution,  we  note  the  following 
geometrical  relationships  for  two  of  the  most  commonly  encountered  head 
closure  designs;  viz.  ,  (i)  full  hemispherical  head  closure;  and  (ii)  full 
ellipsoidal  head  closure. 

(i)  Full  Hemispherical  Head  Closure 


v_ 
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Let  a  and  b  represent  the  semi-major  axis  and  semi-minor  axis  of  the 
ellipsoidal  head  cloture.  Then  the  following  geometrical  relationships  exist: 


Sj2 

co.  v'j 


L  “ 


fcV.bW'' 

— 


R  =  a 

*°  =  90° 

sin  6=1 
ro 

cot  6=0 
ro 

d£  =  R£df 

For  the  case  in  which  the  shell  of  surface  of  revolution  is  a  right  circular 
cone,  the  formulas  for  stresses  may  be  more  appropriately  represented  by 
the  following  expressions. 


3.  2.  c  Uniform  Radial  Shear  Q.  Lb  Per  Linear  Inch  of  Circumference  at 
1,1  —  o  ■■  . 


Q  S%sU<*  (!  2 V  b«r  |t) 

where  C  =  — - - - - - - - - 

We^J.be;  ^  ?.)  +  ^[(be^f.)3 +(be‘i*«j] 

_  Gl0S.S'%*\oA  be?, 5.4  b«.*aS.) 

*>«-**.*  be.'^be^fjt  2v[(b€r,S0)+  (be.'J.)*  ] 

>1  = 
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3.  2.  d  Uniform  Radial  Moment  M  In- Lb  Per  Linear  Inch  of  Circumference 
1  11  —  0  ■■■  — . - . — 


V 

same  expressions  as  for  case  3.  Z.  c  except 

S0 

h  {  *.( 5 . ♦  *  >«•*»  0+2*' [(te^f  .)%  (h*.  ] 


S.ber^  5. 

*****.♦  b€;,^b<O+**[0*nO,v  (b<:t  *.)*]} 
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3.  2.  c  Bending  Stresses  in  Cone  Loaded  by  Uniform  Internal  Pressure  p 

When  there  are  no  edge  shears  or  moments,  but  the  edge  support 
condition  is  such  that  it  is  not  the  same  as  that  necessary  for  the  membrane 
condition,  bending  stresses  will  be  introduced.  The  magnitudes  of  these 
stresses  may  be  obtained  from  the  following  equations: 


ri “  ^ fC  I  [  J  1  *  *»»  )]  + 

c4 b,;^  *  «*".*)]] * 

ro~su'  t  -^r(2b«;it+v5b.:i'5)]  + 

ci ?  5  *  5  (* r 


where 


^>S%-Vo^oirZS.s;«et( 

c  * - - - _ _ 

4-j  S,Cko€v*>I.b«vi5.^  V>eix5.b«^5.)  ♦  2*;[(berJ.)% 

( S.W.' f  ^ zu  be,*^%)  +  b«^5.1 

4{ V»«r^f%t  b« S  s. w.  [(be^ljV  (bo*!.)* ]  j 
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Circular  Plate  with  Uniform  Edge  Moment  In-Lb  Per  Unit  Inch  of 
Circumference 


At  Center: 


At  Edge: 


6(1  -  v)  a  M£ 
*irux ''  eJ7 


9 

o 


12(1  -  v)  aM 


Eh 


T 
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4.  SOLUTIONS  FOR  DISCONTINUITY  SHEAR  FORCES 
AND  BENDING  MOMENTS 


A  tabulation  of  solutions  for  discontinuity  shear  forces  and  bending 
moments  for  pressure  vessel  configurations  of  practical  interest  under 
the  action  of  internal  pressure  is  offered  in  this  section.  The  positive 
sense  of  these  discontinuity  forces  is  defined  in  each  of  the  figures 
associated  with  the  particular  pressure  vessel  configuration  investigated. 

The  following  nomenclature  is  used  in  this  section: 

a  =  >~mean  radius  of  curvature  of  the  circular  cylinder,  inches; 

p  =  internal  pressure,  psi; 

E  =  modulus  of  elasticity,  psi; 

v  =  Poisson's  ratio 

Q  =  uniformly  distributed  circumferential  shearing  force  at 
the  junction,  lb  per  inch; 

s  uniformly  distributed  circumferential  bending  moment  at 
the  junction,  in.  -lb  per  inch 

h  =  shell  thickness,  inch; 

_  3(1  -v2) 

'  T 5T 

1 


Fig.  4. 1.  1.  Two  Long  Circular  Cylindrical  Shells  of  Unequal 
Thicknesses  under  Internal  Pressure 
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J 


a. 


Fig.  4.  1. 2.  Forces  and  Moments  at  the  Junction  of  the  Two  Long 
Unequal  Thickness  Circular  Cylindrical  Shells. 


j 


2Big.  ,  fc  r  ')"=5/2 1 ;  (h,>h2) 


_  (c  -  l)(c  -  1) 

p^2  “  (c*  +  l)2  +  2c3^z(c  +  1) 


:  <hi>h2» 


where 


c  =  thickness  ratio  h^/h^ 


The  above  equations  are  presented  in  graphical  forms  on  Pages  59  and  60. 


J 
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4.  2  Two  Long  Cylinders  of  Unequal  Thicknesses  and  Mismatch  of  the 


Middle  Surfaces 


Fig.  4.2.  1.  Two  Long  Circular  Cylindrical  Shells  of  Unequal 
Thicknesses  under  Internal  Pressure  with 
Mismatch  of  the  Middle  Surface 


-  l)(c5/2  4  1)  +damp1c2(c1/2  4  1) 


(c  +  1)  +  2c  '‘(c  4  1) 


2  -  vf  %w  2  .  1  ,  2.  2  0  1/2  ,  .. 

- *-(c  -  l)(c  -  1)  +  —  da  c  (c  +  2c  +1) 

4pf  2  1X1 


(c  +  If  4  2c  '  (c  4  1) 


M  -  =  M  *  +  —  dpa 
o,2  o,  1  2  m 


where 


c  =  thickness  ratio  hj/h^ 
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Graphs  of  Q  /da  6,p  and  2M  ,/da  d  as  a  function  of  h.  /h?  are  presented 
r  o  in  1  o,  i  n>  i  *> 

on  Pages  61  and  62.  The  Q  and  M_  ,  thus  computed  may  be  superposed  to 

°  O  Op  1 

those  calculated  for  Case  4.1, 


Fig.  4.  3.  1.  Long  Circular  Cylinder  with  Circumferential 
Ring  Stiffeners  under  Internal  Pressure 


J 


Fig.  4.  3.  2.  Forces  and  Moments  at  the  Junction  of  the  Circular 
Cylindrical  Shell  and  the  Circumferential  Ring 
Stiffener 


j 
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M 


Q  =  26M 
o 


4.  4  Long  Circular  Cylindrical  Shell  with  Many  Equidistant  Circumferential 
Ring  Stiffeners 


Fig.  4.4.  1.  Long  Circular  Cylindrical  Shell  with  Many  Equidistant 
Ring  Stiffeners  under  Internal  Pressure 


Q0  =  2pn1BM 


o 


t  / 


-39- 


where 


sinh 

-  sin  6f 

13  : 

’  sinfT 

tr 

+  sin  01 

_  cosh 

b* 

+  cos  (51 

14  : 

’  sinh 

JT 

+  sin  pi 

cosh 

B  t 

-  cos  01 

15  : 

sinh  pi 

-  sin 

Graphs  of  ^4  and  are  plotted  on  Pages  75,  76,  and  77. 
4.  5  Long  Circular  Cylindrical  Shell  with  a  Flat  Head  Closure 

0^  Q0 


Fig.  4.5.  1.  Long  Circular  Cylindrical  Shell  with  a  Flat  Head 
Closure  under  Internal  Pressure.  The  Forces 
and  Moments  are  shown  in  the  Positive  Sense 
above . 
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where 


Fig.  4.6.1.  Short  Circular  Cylindrical  Shell  with  Equal  Thickness 
Flat  Head  Closures  under  Internal  Pressure.  The 
Forces  and  Moments  are  shown  in  the  Positive  Sense 
above . 


J 


(a3  “  a4^  *^1  ~  ^4^15^  *  ^3  ^a5^1  3  *  al^ 

io  =  2pa [(a1-a5nlJ)(b2-b5a13)  -  (a2-a6n14)  <Vb4«i5l 


2  (a4'a3)(b2-b5n13)  +  b3(a2~a6ni4) _ 

°  P  *al  ‘a5ni3^  *b2  "  b5ftl  3^  *  *  2‘a6Qi4*  *bl  "b4ni5*  . 


where  a^  through  a^f  bj  through  b^»and  (5  have  been  previously  defined 
(Case  4.5)  and  ^14*  an<*  ^5  aTe  defined  Case  4.4. 


J 
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4 J  4.7  Short  Circular  Cylindrical  Shell  with  Unequal  Thickness  Flat  Head 

Closures 


Fig.  4.7.  1.  Short  Circular  Cylindrical  Shell  with  Unequal  Thickness 
Flat  Head  Closures  under  Internal  Pressure.  The 
Forces  and  Moments  are  shown  in  the  Positive  Sense 
above . 
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where 


J 


11  h 


(1  -  v) 


1 


A12  =  2(1  -  v) 


A 


-  3a 
13  1 6h  j 


(1  -  v) 


(1  -  v) 


A15  =  2(1  -  v) 


A 


- 

16  ”  16h3 


(l  -  v) 


17 


6(1  -  v) 
02H1h2 


18 


3(1  -  v) 

20Zah2 


_  3(1-  v) 

19 


A 


A 


A 


6(1  -  v) 


20  * 


21  " 


22  = 


_  3(1  -  v) 

20^ah2 

3(1  -  v) 

16p2h2h3 


J 
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B, ,  =  »inhgf  cosh  gl„-  yingf.cpgM  and  .#  ploWed  p#ge  7g 
lb  stall  pi  -  tin  pi 


n17  =  ,inh  fit  .co,.h  Pi  i  »PjL-c-£»-Pl  and  i.  plotted  on  Page  79. 
U  sinn  pi  -  sin4  fit 


„  cot  p 1  sinh  pi  -  sin  pi  cosh  pi  .  .  .  ..  .  _  an 

0  = - 1 - -Xl - X - 1—  and  is  plotted  on  Page  80. 

10  sinh4  pi  -  sin  pi 


n,  q  =  ain  -Pf  ■  £L  and  i8  plotted  on  Page  81 

sinh  pi  -  sin  pi 


2  .  2 

=  g*n^.  JM.  t  sfo-  PJ  and  is  plotted  on  Page  82. 
4U  sinh4  pi  -  sin4  0l 


=  JSllL .  -IP*!  P* - and  is  plotted  on  Page  83. 

sinh  pi  -  sin  pi 


Fig.  4. 8. 1 .  Long  Circular  Cylindrical  Shell  with  an  Ellipsoidal  Head 
Closure  under  Internal  Pressure.  The  Forces  and 
Moments  are  shown  in  the  Positive  Sense  above. 
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The  above  expressions  for  Qq  and  are  plotted  as  functions 
on  Pages  59  and  60. 
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Fig.  4.9.  1.  Short  Circular  Cylindrical  Shell  with  Equal  Thickness 
Ellipsoidal  Head  Closures  under  Internal  Pressure. 
The  Forces  and  Moments  are  shown  in  the  Positive 
Sense  above. 


•47. 


where 
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4.  10  Short  Circular  Cylindrical  Shell  with  Unequal  Thickness  Ellipsoidal 
Head  Closures 


M 

M 

Q 

Q 


1 


2 


l 


2 


Fig.  4,  10. 1*  Short  Circular  Cylindrical  Shell  with  Unequal 
Thickness  Ellipsoidal  Head  Closures  under 
Internal  Pressure.  The  Forces  and  Moments 
are  shown  in  the  Positive  Sense  above. 
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A> =  - 


ai2  ’  P-JB7  V31' -  “  1 


The  functions  ^  through  ft  are  identically  defined  as  those  for  Case  4.7 


4.11  Long  Circular  Cylindrical  Shell  with  a  Hemispherical  Head  Clo»ure 

ft.  GU 


Fig.  4.11.1.  Long  Circular  Cylindrical  Shell  with  a  Hemispherical 
Head  Closure  under  Internal  Pressure.  The  forces 
and  moments  are  shown  in  the  positive  sense  above. 
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c  c 


ZEVs 
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where 


3Q-y*) 
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4.12  Short  Circular  Cylindrical  Shell  with  Equal  Thickness  Hemispherical 

Head  Closures 


a.  a,  q. 


a*  Q.  a.  a. 


Fig.  4.  12.  1.  Short  Circular  Cylindrical  Shell  with  Equal 
Thickness  Hemispherical  Head  Closures 
under  Internal  Pressure.  The  Forces  and 
Moments  are  shown  in  the  Positive  Sense 
above . 

The  formulas  for  Q  and  M  for  this  case  are  identical  to  those  for  Case  4.9. 
o  o 

The  dimensionless  influence  numbers  a,,  a.,,  •  •  •  ,  blt  *  •  •  •  b-  etc., 

1  C  1  D 

are  also  identical  with  the  exception  of  a.y  where  for  this  case, 

a3s-?°  "V) 


Fig,  4. 13. 1.  Short  Circular  Cylindrical  Shell  with  Unequal 
Thicknet i  Hemispherical  Head  Closures 
under  Internal  Pressure.  The  Forces  and 
Moments  are  shown  in  the  Positive  Sense 
above . 

The  formulas  for  Q^,  and  for  this  case  are  identical  to  those 
for  Case  4. 10.  The  dimensionless  influence  numbers  and  must, 
however,  be  changed  to  the  following 


Fig.  4. 14 . 1.  Long  Circular  Cylindrical  Shell  with  a  Conical  Head 
Closure  under  Internal  Pressure.  The  Forces  and 
Moments  are  shown  in  the  Positive  Sense  above. 


j 
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Graphs  of  bj,  and  a s  a  function  of  and  o  for  v  r>f  0.  30 

arc  plotted  and  arc  presented  on  Pages  &6throtigh  #9. 

4.15  Short  Circular  Cylindrical  Shell  with  Unequal  Thickness  Head  Closures 
of  Ellipsoidal  and  Conical  Shape 


Fig.  4.  15.  1.  Short  Circular  Cylindrical  Shell  with  Unequal  Thickness 
Head  Closures  of  Ellipsoidal  and  Conical  Shape  under 
Internal  Pressure.  The  Forces  and  Moments  are 
shown  in  the  Positive  Sense  above. 
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The  functions  0^  through  arc  identical  to  those  defined  in  Case  4.7. 

4.  16  Short  Circular  Cylindrical  Shell  with  Unequal  Thickness  Head  Closures 
of  Hemispherical  and  Conical  Shape 


Fig.  4.  16.  1.  Short  Circular  Cylindrical  Shell  with  Unequal  Thickness 
Head  Closures  of  Hemispherical  and  Conical  Shape 
under  Internal  Pressure.  The  Forces  and  Moments  are 
shown  in  the  Positive  Sense  above. 


The  formulas  for  Q^f  Mj,  and  for  this  configuration  are  identical  to 
those  for  Case  4.  15  except  that  is  now  defined  by 

AI3  = 
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(it-Z)  d 
*t )'tfZ 


Fig.  5.2  Variation  of  Juncture  Shearing  Forces  wi*  T^'ckness  o  for 
Two  Long  Cylindrical  Shells  of  Uneq  ^1  Tv  ckness 
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Fig.  S.  3  Variation  of  Juncture  Shearing  Forcee  with  Thicknesft  Ratio  Due  to 
Mismatch  of  Middle  Surfaces  Only  of  Two  Long  Cylinders  of  Unequal 
Thickness  (to  be  used  in  conjunction  with  Fig.  5. 1} 
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Fig.  5.5  Variation  of  ft.  with  f3x 
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Fig.  5. 8  Variation  of  with 
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Fig.  5.9  Variation  of  with 
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Fig.  5.13  Variation  of  with  pi 


71 


COSH 


0  2.0  4.0  6.0  8.0  10.0 

A* 


Fig.  5. 14  Variation  of  with  pi 
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Fig,  5,15  Variation  of  fl^  with  |3i 
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Fig.  5.17  Variation  of  R.jWithf)i 
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Fig.  5. 18  Variation  of  with  pi 


Fig.  5. 19  Variation  of  with  §/ 
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Fig.  5.21  Variation  of  o17  with  p| 
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Fig.  5.  22  Variation  of  fyg  with  |5f 
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Fig,  5,  24  Variation  of  n2Q  with 
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Fig.  5. 26  Variation  of  n22  with 
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Fig.  5.  27  Variation  of  ^3  with  pi 
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Fig.  5. 28  . Variation  of  cosec  a  with  6 


Fig.  5.29  Variation  of  cot  a  coaec  a  with  s 
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Fig.  5. 30  Variation  of  4b ^  cot  a  with  4 
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APPENDIX 


DERIVATION  OF  FORMULAS 


A.  BENDING  STRESSES  IN' SHELLS 

1.  Symmetrical  Deformation  of  Circular  Cylindrical  Shells 

The  differential  equation  governing  the  symmetrical  deformation 
of  circular  cylindrical  shells  in  accordance  with  the  classical  shell  theory  is 


(A. 1.1) 


where  D  =  flexural  rigidity  of  the  shell 


12(1  -  v2) 


Z  =  load  intensity 
h  =  thickness  of  shell 
a  =  mean  radius  of  curvature 
E  *  Young’s  modulus  of  elasticity 
v  =  Poisson*  s  ratio 
ur  r  radial  displacement 

In  the  case  of  uniform  shell  thickness  Eq,  (A.  1.1)  becomes 


(A.1.2) 


4u 

2  r 
a 


=  Z 


A-l 


(A.  1.3) 


Using  the  notation 


p4  =  &  .  iiijja 

4a2 D  a2h2 


(A.  1.4) 


Equation  (A.  1. 3)  can  be  represented  in  the  simplified  form 


d4u 


T  +  4P4«r=^ 


dx 


(A.  1.5) 


for  which  the  solution  is 

ur  =  e^X(C^  c os  px  +  sin  Px)  +  e  ^X(C^  cos  Px  +  sin  Px)  +  f(x)  (A.  1.6) 

In  Eq.  (A.  1.6),  f(x)  is  a  particular  solution  and  C^,  C^»  C^,  and  are  the 
constants  of  integration  which  must  be  determined  in  each  particular  case 
from  the  conditions  at  the  ends  of  the  cylinder. 

The  expressions  for  the  membrane  forces,  bending  moments,  and  transverse 
shearing  force  associated  with  the  symmetrical  deformation  of  circular 
cylindrical  shells  are 


(cont.) 


A-2 


M9  =  -  VD— 
W  dx 


Qx=-D- 


(A.  1.7) 


In  the  case  where  there  it  no  pressure  Z  distributed  over  the  surface  of  the 
shell  and  if  the  end  condition  is  such  that  Nx  =  0,  then  f(x)  =  0  in  Eq.  (A.  1.6). 
Furthermore,  for  long  cylindrical  shells  subjected  to  loading  conditions  such 
that  for  large  positive  values  of  x,  the  deflection  is  finite,  :  0  and 

Eq.  (A.  1.6)  reduces  to 


u  =  e’^IC,  cos  px  +  CA  sin  Px) 
r  3  •k 


(A.  1.8) 


For  Nx  =  0,  the  first  expression  of  Eq.  (A.  1.7)  gives 


du  u 

x  _  v_r 

ax  a 


which  when  substituted  into  the  expression  for  yields 


Ne=-  — 


(A.  1.9) 


The  stresses  in  the  shell  are  given  by  the  expressions 


<r  =  ± 
x 


Nfl  6M. 

vtV 


(A.  1.10) 


where  <r  =  meridional  stress 

<Tq  =  circumferential  stress 

2.  Coordinate  System  and  Signs  Convention 

The  following  signs  convention  for  the  rotation  and  deflection  are 
used  throughout  this  report: 

Positive  deflection  -  radially  inward  with  respect  to  the  center- 
line  of  the  cylinder 

Positive  rotation  -  clockwise  viewing  the  upper  cut  of  the 
cylinder. 

The  signs  convention  for  the  forces  and  moments  are: 

Positive  shearing  force  -  upward  on  the  left  end  and  downward  on 

the  right  end  of  the  shell  element,  the 
left  end  being  at  the  plane  of  the  origin. 

Positive  moment  -  tensile  on  the  inner  surface  or  com¬ 

pressive  on  the  outer  surface  with 
respect  to  the  center -line  of  the  cylinder. 

The  above  signs  are  graphically  represented  in  Fig,  A.  2.  1. 


J 


A-4 


i 


In  the  case  of  a  long  circular  cylindrical  shell  under  the  action  of  uniform 
radial  shear  Qq  lb/ in*  along  the  circumference  at  the  end,  the  two  constants 
of  integration  and  in  Eq.  (A.  1.8)  are  determined  from  the  conditions 
at  the  loaded  end;  viz.  , 


Specifically,  they  are 


(A.  3,1) 


(A. 3.2) 


(A. 3. 3) 


These  values  of  and  may  now  be  substituted  into  Eq.  (A.  1. 8)  to  give 


u 

r 


Q  e  cos  Px 
o 

203D 


The  successive  differentiations  of  Eq.  (A.  3.4)  are 


2PZD 


(sin  px  +  cos  px) 


(A. 3. 4) 


(cont  ) 


A-6 


d*u  Q  e’^x  sin  px 
r  o 


d3»  o  .-*■ 

X  O 


(cos  Px  -  sin  Px) 


(A.  3.  5) 


The  slope  is  therefore  given  by 


du  Q  ft 

—t~  =  — x~~  e”"x  (sin  px  +  cos  Px) 
dx  ZpD 


(A.  3. 6) 


and  at  the  loaded  end,  x  =  0;  therefore, 


2P2D 


(A.  3.7) 


The  expressions  for  the  bending  moments  M^,  Mg,  and  the  transverse  shear 
force  Qx  are,  in  accordance  with  Eqs.  (A.  1. 7),  respectively, 


U  ft— 

Mx  =  -jgp  c  sin  Px 


vu  a 

w  o  -px  .  a 

Mg  =  — j p—  e  r  sin  Px 


Q  =  O  K  fcoa  0x  -  ain  flx) 
x  o 


(A.  3.  8) 


and  the  hoop  force,  Nfl,  according  to  Eqs.  (A.  1.9),  (A.  3.4),  and  (A.  1.4) 


becomes 


Na  =  2Q  Pae*Px  cos  Px 

v  O 


(A.  3.  9) 


The  expressions  for  the  stresses  in  the  ehell  aie  therefore 


r  <inPx 

x  Ph2 


<r#  =  -j—  e*^*  (Pa  co«  Px  ±  ein  Px) 


(A.  3. 10) 


The  maximum  deflection  ie  at  the  loaded  end,  where 


Ufmax  =  Ur(x=0) 


(A.  3. 11) 


Let  us  now  define 


flj  «  e  sin  px 


(A.  3.12) 


Rj  =  e"^x  cos  Px 


(A.  3. 13) 


Rj  =  e*^x  (sin  Px  +  cos  Px) 


(A.  3. 14) 


Then,  Eqs.  (A.  3. 10)  become 


<r  =  ± — k  fl, 
x  Ph1  1 


_ o 

*9s-r 


(p*^  *  w  °i) 


(A.  3. 15) 


J 


and  Eqs.  (A.  344),  (A.  3.  6)  will  become,  respectively, 


a  = 
r 


(A. 3. 16) 


(A,  3,  17) 


In  Eq.  (A.  3. 15)  where  ±  signs  occur,  the  upper  sign  refers  to  the  stress 
state  at  the  inner  surface  of  the  cylinder  and  the  lower  sign  to  the  stress 
state  at  the  outer  surface  of  the  cylinder. 

4.  Long  Circular  Cylindrical  Shell  under  the  Action  of  Uniform 
Radial  Moment  MQ  in.  -lb/ in.  of  Circumference  at  End 


In  this  case,  the  constants  and  in  Eq.  (A.  1.8)  are  determined  from  the 
following  boundary  conditions: 


Upon  application  of  these  boundary  conditions  to  Eq.  (A.  1.8),  there  results 


M 

o 

2P2D 


(A. 4. 2) 


With  these  values  for  and  C^,  the  expression  for  the  deflection  Ur  becomes 


u 

r 


m  « 

— c  PX  (sin  px  -  cos  px) 
ZpD 


(A.4.3) 


The  successive  differentiations  of  Eq.  (A.4.3)  are 


du  M 
_ _ r  o 

dx  =FD 


cos  Px 


(sin  Px  +  cos  Px) 


d3u  2M  » 

- y  =  — jy  Pe  sin  Px 

dx 


(A.  4.  4) 


.giving  for  the  slope 


9  = 


du  M 
r  o 

“ST  =  ? TT 


■px 


COS  Px 


(A. 4. 5) 
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A-10 


The  slope  at  the  loaded  end  where  x  *  0  is.  therefore. 


M_ 


e(x=o) s  pur 


(A.  4.6) 


Equations  (A. 4.  3)  and  (A.4.4)  may  now  be  introduced  into  Eqs.  (A«  1. 7)  and 
(A.  1. 9)  to  give  the  following  expressions  for  Mx<  Mg,  Qx  and  Ng.  Indeed, 


M  =  M  e  (sin  Px  +  cos  (be) 
x  o 


Mft  =  vM  e"^x  (sin  Px  +  cos  Px) 

U  O 


O  =  -2M  pe  *^X  sin  Px 
x  o 


Ng  =  2MQP^ae  ^  (cos  Px  -  sin  Px) 


The  expressions  for  the  stresses  in  the  shell  are,  accordingly, 


6M 


<r  =  *— r*  e  ^X  (sin  px  +  cos  Px) 


2M 


s-Px  |p2a  jco#  px  _  sin  pxj  ±  3v_  px  +  cos  px^J 


or,  by  introducing  the  notations  of  Eqs.  (A.  3.  13)  and  (A.  3. 14),  vix. 


=  e~^x  cos  Px 


(A.  3. 13) 


flj  =  e  ^X  (sin  Px  +  cos  px) 


(A.  3. 14) 


A- 11 


and  defining 


^4  =  e  (co»  Px  "  **n  P*) 


we  find 


and 


6M 

h 


2M  ,  - 

U  ~TT (P  an4  *TT  n3) 


(A. 4. 7) 


(A. 4. 8) 


(A. 4. 9) 

(A. 4. 10) 


(A. 4. 11) 


5. 


Short  Circular  Cylindrical  Shell  Bent  by  Force,  Distributed 
Along  the  Edge, 

(a)  Equal  Edge  Force*. 

<*o  Q. 


In  the  case  of  shorter  shells  loaded  such  that  no  pressure  Z  is  distributed 
over  the  surface  of  the  shell,  the  deflection  equation  (A.  1.  6)  can  be  put  into 
the  following  form  by  the  introduction  of  hyperbolic  functions  in  place  of  the 
exponential  functions.  Thus, 

Tly.  ®  ^  C|Ci^pXsUV\|ix  +  C^Cu^XCi^pX 

(A.  5.1) 

By  selecting  the  origin  of  coordinates  to  be  at  the  middle  of  the  cylinder,  it  is 
readily  seen  that  Eq.  (A.  5*  1)  must  be  an  even  function  of  x.  Accordingly, 

Cz  -  C3  =  0  (A.  5. 2) 


and  Eq.  (A.  5. 1)  reduces  to 

=  C|  sin  fix  sinh  fix  +  cos  (he  cosh  Px  (A,  5.3) 


A-13 


A.14 


The  successive  differentiations  of  Eq.  (A.  5. 6)  are 


dfuT  m  2Q, 

dx*  “  '  O 


(A. 5. 7) 


The  deflection  and  the  slope  at  the  loaded  ends  are  therefore,  respectively, 


and  the  expressions  for  M^,  Q^,  and  become,  respectively. 


20.,  c*s^t0wW§s:'*f  x.s  «**Upx  •  sU  v4 

^  Sin  ^  ♦  JinUp^ 


(cont.) 
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2l>(X,  C#*  mW^  si  Il-Wpx  —  S‘<n^si«V,||>e»*^XC*lU^S 
P  S{«jU+*S«W^ 


(A.  5. 9) 


61  -  -20. 

x  • 


Cos-^  <«W^  (c.l  fxtinUfx  4  Sm^XC^px)-  Sln^sUi^f  (on^X«i»Up*  -Si't^X  OmU^ 

Sin^l  +  S>»Vi^J 


40#pa  (*4f.aJ£ 

Stv^pX  t  S%nV»  p€ 


Hence  the  stresses  in  the  shell  are,  respectively, 


(T  »  T- 

x 


St  w  *' -  9 '«v§£  Cr»px 0+^pX  j  ^  ^ 

- I  <A* 

^  •  Sin^t  "fr  *»nkp  J 


v- 


^  Q 

— ‘ - ! - -  ±  TT  (CoS§  e*‘u^  *:"P*  *'~u?x  ‘ •**•? j;-u^2  c«p*o«kpx) 

^in^/  +  s!*V(^J  P 

“  ^s'"^  4  Cd^xto 


lat  us  now  define  the  following  parameters: 

=  sin  Px  sinh  Px 
n6  =  sin  px  cosh  px 
ilj  =  cos  Px  sinh  px 


(cont.) 
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fig  =  co*  Px  co»h  px 


Pi  .  .  Pi 

•in  -jr  *inh  -jr 

^9  =  .TnUTTOT 

co*  co*h  ^ 

ni2 s  TnrprrjiiSiw 

-  »inh  Pi  -  tin  pi 
n13  “  sin  pi  f  *inh  pi 


co»h 
n14  ~  Tin 


i  Pi  4  cos  Pi 

p?  r.iafr 


(A.  5.11a) 


(A.  5.11b) 


Then  the  expressions  for  the  stresses  in  the  shell  become,  respectively, 


1ZQ 

T  <°8«9  '  W 


{A.  5. 12) 


«rQ  =  [  PatRgl^  +  Cgf^)  ±  (A.  5. 13) 


and  the  expressions  for  deflection  and  slope  are,  respectively, 


aQ 

«¥»»  +  W 


(A. 5. 14) 


(A.  5. 15) 


which  become,  at  the  loaded  ends. 


u  • 

r(x=*|) 


(A.  5. 16) 


A-i7 


(A.  5.  17) 


In  this  case,  the  four  constants  of  integration  of  Eq.  (A.  5.  1)  will  be  evaluated 
from  the  following  boundary  conditions: 


J 


A-18 


Specifically,  they  are 


C,  =  O 

_  (g^Cwpl-  <aa»;»upj) 

1  (sCnUy  - 

(ati;»up*  -  at*~frQ 

3  2p30  (zUl? -  SJ«1^) 

£  _  Q|  (si>>jS^  ~*  Co»V^t^4 

4  2  £3D  (shvW^jf  -  tUty!) 

(A. 5. 19) 


These  values  of  c .  (i  =  1 , 2,  3,  4)  may  now  be  substituted  into  Eq.  (A.  5.  1) 
to  give 

_ _ i _ 

r  2  p3E>  (*UV>*|Si  -  *:«*pX) 

(Q1s;«p{-  a,3U,pt)  s'-v^U  o^^x»;*upx.  + 

1 

[<Jz(c*sp£»;^f£-  *up£e%»V,^)f  fll^suVft  or,U|lt-*;^<<>«(U)]o-pxc«t»p*. 


(A.5.Z0) 


^QgSmUpf  -  Q ,  sUflt  ^  *» ■*  p*  C*tV»P*.  + 


A-19 


The  successive  differentiations  of  Eq.  (A.  5.  20)  are 


<K  . 

obt 


.*  ~f  1 ,  - ”  |  14  c«*^e-V^x) 

2|3  D(^«U  fX-sUdl)  I 

+  s:^pi  (at  supi  - »^u^)(c^3c  c*v^*.  -  »;^xj'-xV,p>c) + 

(  c^|S**i-V>px  -  *n*'f*)^*(e*tf  *  »UV»jU  -  tU|J(  )  + 

Q(  (s'^UfJ t  -*in  ft  ®*»/**)J  | 


(A.  5. 21) 


dx* 


(afcs»*.»W^{  -  <?(  *U  jU  )  c.*px  slhVfx  — 
suufi  al*'.nU^t)»i«.p»o»»vpx  - 


PDfrjffJl-tU'ft) 


d*3 


- ! - 1  s.~gt 

—  0f»Upxj|l9i^Cu^t StitK^  ^ 

Q^-wpt  o*>kpt-sU 


The  deflection  and  the  slope  at  the  loaded  end  where  x  s  0  are  then 


(A. 5. 23) 


Q  —  -  +  s:«v,V) 

z (suUTpjf- 


On  the  other  hand,  the  deflection  and  the  elope  at  the  loaded  end  where 
x  =  f  will  be 


■^(x-l)  =  -  Q_»(V,,>V^ “»»?*) 4  O, 

■2  {?& (s ;  -  si*^) 


(x»  ( )  z  '  "  - - 

2f  D(s;KV,lfe-sU^Jt) 

The  expressions  for  M^,  M#,  Q^.  and  N#  may  now  be  easily  obtained  by 
•ubetituting  the  expressions  for  uf.  d2ur/dx2,  and  d3u  /dx3  from 
Eqs.  (A.  5.  20)  and  (A.  5.  21)  into  Eq.  (A.  1.  7).  Of  particular  interest  to 
our  study  are  the  stresses  in  the  shell.  These  are  found  to  be 


(A.  5. 24) 

(A.  5. 25) 
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1  K(*u.*fU- *••»>*  ) 


J 


+  '*^p»e«u^t) 

+  QUc»»p*c.*Wj}xT”  *'^^xa;-V,px)J^(Jl(o»»pt3.s^-*:-.^i  + 


Q((sUU^<  CwUpt  -*:«^£  c»«|?Jt)J  ’ 


Using  the  parameters  tt$,  0^,  tty,  fig  defined  in  Eq.  (A.  5.  11)  and  introducing 
the  following  additional  parameters: 


A 


IC 


A 


H 


20 


A« 


A 


2Z 


XI 


*5 


s1hV,'(U  -  s:*Y«. 

Crtjlt  sUVi^t  ~  S'np(  Cwti^t 

suv>ap*  - 

S»nU  jlf-  *"  3»n  P't 

s;^(^  -  *^yi 
*uu'M-  »u*f* 

-  s**xpZ 


(A. 5.26) 


J 
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We  find  the  expressions  for  the  streseee  in  the  shell  to  be 


(\vv -)AS] 

°i  “  t|?a[A«(a^al-A„sl)  +  il7(A,lat-i^lQt)+ 

*  *jjv7  [-^7  O^wV 

■^“6 1  ® l)  "  5 CA»  'A  14  ®  I  )J^  (A.  5. 28) 

The  expressions  for  deflection  and  slope  are  correspondingly 

(-^1*^1* -&•*<*  J^g]  (A.  5. 29) 

®  iX -^S^b)  +  (-^,ri<az~Az^,)(ilg-{l5) 

+  (A.  5.  30) 

and  at  the  loaded  ends,  the  deflections  and  slopes  become 

^(a-o)-  “  (A.  5. 31) 
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■^(*•0“  <A-5-3Z) 

(*  ‘^*‘*1®*”  |)  (A.  5.  33) 

Vo  "  "  7f7 (il’ ®1’  2Jl",a,)  ,A'5-34' 


6.  Short  Circular  Cylindrical  Shell  Bent  by  Moments  Distributed 
Along  the  Edges 

(a)  Equal  Edge  Moments 


Again,  in  this  case,  Eq.  (A.  5.  3)  is  applicable  The  constants  and  c^ 
may  be  determined  from  the  following  boundary  conditions: 


v.-i)-  -o(^4=  o 


j 
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Specifically,  they  are 


M#  /Co*4£  s*»«v4| 

*  0*0  V  *uvpi4-»up t  / 

c  M.  /C*^  3^4  ~  \ 

^  j5lt>  \  »i«U^4*«*npt  ) 


(A.6.Z) 


The  deflection  equation  for  the  middle  surface  of  the  shell  is  then 


M0  4 « ' «J c  stp^p*.  f  (c^ 

^  D  41 


(A. 6, 3) 


The  successive  differentiations  of  Eq.  (A.  6.  3)  are 


2M#(c«s^s>mU^  c«^»^*px+st.^£c9iU^ 


E>(sl*Vft+»^pi) 


A  ^MJ^aU^cwU^cwpxtuU^x- 


<A.  6. 4) 
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The  deflection  and  slope  at  the  loaded  ends  are  now  obtained  by  substituting 
±  t/Z  for  x  in  the  above  equations.  Specifically,  they  are 


-j 


0,  _ 

(*“  *  *)  f  ©  \  5,'^v  04  4  / 


|U  4  *:«|U 


(A.  6. 6) 


The  stresses  in  the  shell  are,  respectively. 


5J~U  * + *-4 

4  -  s;^c^^^^ac%tp^©^px  4^  v*^xs:Jkpx)| 


Using  the  parameters  12,.,  12^,  flg,  defined  in  Eq.  (A.  5.  11),  and 

introducing  the  following  additional  parameters: 


$ 


IL 


!* 


tosU-^ 

siv>  pX  ^  pX 


Cm^ 
Sm|U  4 


(cont.) 
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(A. 6. 7) 


SL 


(tf  4 


We  find  the  streasee  in  the  ehell  to  be 

A  (A.6.8) 

^  vT  ^  ^*e)  +  ^8)J  * 

(A.  6.9) 

together  with  the  following  expressions  for  deflection  and  slope: 

-•v -  -  [( Aj  a5*  (a,,-  a.)  a,]  <*•  *■ » 

9  “  •A-6  n: 


The  corresponding  expressions*for  deflection  and  slope  at  the  loaded  ends 
are  then,  respectively. 


~U 


(*•*■£) = 


m.h,. 


(A.  6. 12) 
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V-*4)  ~  *  ‘A-6-‘3' 

(b)  Unequal  Edge  Momenta 


M 


l 


Here,  because  of  the  unsymmetrical  loading,  Eq.  (A.  5.  1)  must  be  used, 
and  the  constants  of  integration  will  be  determined  from  the  following 
boundary  conditions: 


J 
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Specifically,  they  are 


To  recapitulate,  the  deflection  Eq.  (A.  5. 1)  is 


VT  —  4 


with  eucceesive  derivatives  of 

“  f*[C<  4  C^eSKjbtfj^px+CwpxewU^x)-!- 

C^Cc^pxft-kpx-  4  C4(eoipx*i-kpx- «:,p*owVp*)J 


jO  -  xcoV^t  ♦  Ctc*ipx*J-.V^t-  Cj*I->p*.o«*'P>t-C^»u^cfU^,px) 


1  ~  “  Z^fc^Co^xli^X*  Sv^XCtA^x)  4  C^CM^XOnl^H-tti^XSiJ^X^' 
(Ax  L  «k 

Cj(s:>fxsi.k^l  VCrt^xCuV^x).  C^(u*ptOnLpX  4  CtfU «lJ^x)J 
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The  deflection  end  the  a  lope  at  the  loaded  end  where  x  =  0  are  then 


J 


(x.»  a) 


-  C, 


.  (A.6.l6) 

2  (slnu'ft-  S 


_  KA ,  (s;^fc ♦  suvftf  c  oV>pl)  -  VAz(*'*ft  +  ) 

|&D  sU’ptf) 


(A.  6. 17) 


and  the  deflection  and  elope  at  the  loaded  end  where  x  =  I  are 


^ f*  s  i*Mt  ♦  Ct(s^pf  oW^f  ♦  e  .»fX  »Ul,^£)  ♦  C+coaft  CvXft 

«  >  f  2M,  »W<  "j 

zpzo  I 


J 


«uuxpe  -  s.v'pA 


j 


(A.  6. 18) 


0(ot*£)  ®  ^[C,(^|0e^p|  +  C**p*5U,f*)i-  ZCt©«*p|c*»),p£ 

SU  V>p ^ 


p 


\  M|  (Smpto»»V,ff +  0^*  sUU^)  -M2(jUpt  Cffpt  4s.n\.ffto«A,^)' 

s^w'pi-ss^at 


(A.  6. 19) 


J 
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The  stresses  in  the  shell  are*  respectively, 


,  co»^x  e«*|hc  *  Clo»p-x SUU^X  -  C,  s;^x  coUffcx  - 

m 

6  C«Vp(  t  CA^-  1 

*\L  suw,lfi  -  »UzJl<  J 

C^VjSx-  OapxaUWpx!)  V  M,  C*i|?x  toUjSx  — 

L  J  r  r  |A-‘-J0> 


*» 

2  fa  M  (»U.*| it ♦sinSt)  -  2NLs5~M  SJ4.M 

u  [  y^-*«Kx|te  T  n 

v  Ma(»uftc^i.pC4  c 


Cos^xsUVx^))  t 


iK^XC^Wp 


|A.fr.  21) 
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Using  the  parameters  0^,  0^,  0^,  0g  defined  in  Eq.  (A.  5.  11),  and  82^, 
02  j  defined  in  Eq.  (A.  5.  26)  and  introducing  the  following  additional 
parameters: 

^  4  a«r»(l X 

*7  -  f  i 


&*>Vy  | 

5^V,*pl  -  sUa(U 


(A.  6. 22) 


We  find  the  stresses  in  the  shell  to  be 


±  J*[(Ml%“Mx£lB)CILc-il7)"(MiA*.“2^»,)^B*Mia8J  (A. 6. 23) 


^  (A.  6. 24) 

The  expression  for  the  deflection  and  slope  then  becomes 

Vt  “  Tp^O  1  ^ Q  »X-a*6+'Q’7) "  (M *  ^.~'2 M  t  )  A#-  M ,  A*]  (A.  6. 25) 


-M.CA.4A7)] 


(A.  6. 26) 
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At  the  loaded  ends,  the  expression  for  the  deflection  and  the  slope  becomes 


^ (A.  6.27) 
^<^-1)-  (A.  6. 28) 

""  (A.6.29) 

V"0  ~  (M.-^o“M»^»7)  (A.6.30) 
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B.  JUNCTURE  SHEAR  FORCES  AND  BENDING  MOMENTS 


Fig.  B.  1.1.  Two  Long  Circular  Cylindrical  Shells  of 

Unequal  Thickness  Under  Internal  Pressure 

Under  the  action  of  internal  pressure,  the  circular  cylindrical  shells  will 
extend  radially,  each  by  a  different  amount  due  to  the  different  thicknesses. 
However,  in  the  actual  vessel,  the  two  cylinders  are  kept  together  and  strain 
compatibility  at  the  junction  indicates  that  there  must  act  shearing  forces  Qq 
and  bending  moments  Mo  uniformly  distributed  along  the  circumference  and 
of  such  magnitudes  as  to  eliminate  this  continuity. 


A-34 
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Eh. 

D= - L—  ;  (i=1.2) 

1  12(1  -v6) 

a  =  mean  radius  of  curvature  of  the  circular  cylinders; 
p  =  internal  pressure,  psi; 

E  =  modulus  of  elasticity,  psi; 
v  =  Poisson's  ratio 


(B.  1.5) 


Qq  =  uniformly  distributed  circumferential  shearing  force  at  the 
junction,  lb  per  inch; 

Mq  =  uniformly  distributed  circumferential  bending  moment  at  the 
junction,  lb  per  inch; 


J 


Fig.  B.  1.  3  *  Signs  Convention  for  Middle  Surface  Radial 
Deflection  and  Rotation  at  Junction 

For  strain  compatibility,  we  must  have 

6  (for  Cylinder  No.  1)  =  5  (for  Cylinder  No.  2) 

6*  (for  Cylinder  No.  1)  =  6*  (for  Cylinder  No.  Z)  (B.  1.6) 

J 
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Hence 


(B.  2.4) 


L.. 


s'“  ***s.) 


where  j  and  Mq  2  arc  the  uniformly  distributed  circumferential  bending 
moments  at  the  junction  for  cylinders  1  and  2,  respectively,  and  are  related 
to  each  other  by  the  following  equation: 


M  «  =  M  ,  -  4-dpa 
o,  1  o,  2  2  r  m 


(B.2.  5) 


where 


*m  IU1  ♦  *2> 


(B.  2.  6) 


If  we  now  substitute  Eq.  (B.2.5)  into  Eqs.  (B.2.3)  and  (B.2.4),  we  obtain 
for  cylinder  No.  2, 


(B.2.8) 


Again  for  strain  compatibility,  we  must  have 

6  (for  Cylinder  No.  1)  =  6  (for  Cylinder  No.  2) 
6*  (for  Cylinder  No.  1)  =  6'  (for  Cylinder  No.  2) 
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Hence, 

J 

[v^K^r-o] 

from  which 


a.  = 


M..r 


_2 

ZE 


«L  .  «1Y-L  ♦  _!_)  +  _*a-  f_L  4  ±) 
=  \  k.  w.  fw  f.  pJ 

_l(J_  +  _L-\%.  _1 _ /_L  .  _L\ 

2W>.  tf*V  fl?/ 

w«i  ^vj _ l\  a/j_ + — * —  *  _j_  \ 


f>  (B.2.9) 


f>(B.2. 10) 


Equations  (B.  la.  9)  and  (B.  la.  10)  may  be  put  into  the  following  convenient 
forms: 

■|j-(c-ilXc*Y4 *  I) 

Qo  —  - — - - , - -jp 


J 
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where  c  =  thickness  ratio  hj/h^ 

i  =  middle  surface  radius  ratio  Sj/a^ 

Observe  that  the  above  equations  are  developed  for  aj<  a^.  If  a^  >  a^>  then 
d  becomes  negative,  and  the  components  of  shear  and  bending  moment  due  to 
the  mismatch  alone  of  the  middle  surfaces  change  in  direction. 

For  thin  shells,  a/h  is  large  and  i  can  be  taken  as  unity  for  all  practical 
purposes.  The  equations  for  Q  ,  M  ,  then  simplify  to 

O  Of  1 


^(cmXc*»|)+«U,JS,cVM) 

— Li - - -  ip 

(c*+  l)*  +  2cV‘(c*0 

fg-(c-iX^-0v^(cV2^i) 

(c*+  lj)*+  2cV'(c+  I)  r 
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Fig.  B.3.  1  -  Long  Circular  Cylinder  with  Circumferential  Ring 
Stiffeners  under  Internal  Pressure 

Basic  Assumptions:  (a)  The  ring  stiffener  spacing  is  such  that  the  influence 

of  one  does  not  extend  to  the  next. 

(b)  The  term  "a"  will  be  taken  to  be  the  outside  of  the 
shell  instead  of  mean  radius.  For  shells  with  a 
large  radius  to  thickness  ratio,  the  effect  is  negligible. 


Fig.  B.3.2  -  Forces  and  Moments  at  the  Junction  of  the  Circular 

Cylindrical  Shell  and  the  Circumferential  Ring  Stiffener 
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Boundary  Conditions:  (a)  The  radial  deflections  of  the  shell  and  ring  are 

equal  at  the  junction, 

(b)  The  slope  of  the  shell  adjacent  to  ring  it  zero. 

Analysis:  The  deflection  of  the  long  circular  cylindrical  shell  at  the  junction 
produced  by  the  loads  p,  Qq  and  Mq  is,  from  Reference  4, 


(B.3.1) 


The  radial  displacement  of  the  circumferential  ring  stiffener  under  the  action 
of  the  forces  2Qq  and  p  may  be  readily  obtained  from  the  plane  stress  solution 
of  a  thick-walled  cylinder  under  the  action  of  internal  pressure.  Specifically, 
at  the  ring-shell  interface, 


.,lOo  ♦  pc)  A2  t  2  \ 

8 - ^1777*7 


(B.3.2) 


In  accordance  with  boundary  condition  (a),  we  can  immediately  write 
,  2  a(2Q  +  pc)  /,  2  2  \ 

=  — k — (77 7*j 

The  slope  of  the  circular  cylindrical  shell  adjacent  to  the  ring  produced  by 

the  forces  p,  Q  and  M  is,  from  Reference  4, 
o  o 


6*  =  — 5— (2PM  -  O  ) 

zpV  °  ° 


(B.  3.4) 


This,  in  accordance  with  boundary  condition  (b)  is  zero.  Hence 


Q  =  2PM 
o  o 


(B.3.5) 
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J 


Equation  (B.3.  5)  may  now  be  substituted  directly  into  Eq,  (B.  3,  3)  to  give  an 
explicit  expression  for  the  circumferential  bending  moment  Specifically, 


Q  =  23M  (B.3.7> 

o  o 


Fig.  B.4. 1  -  Long  Circular  Cylindrical  Shell  with  Many  Equidistant 
Ring  Stiffeners  under  Internal  Pressure 

Basic  Assumptions:  (a)  The  ring  stiffener  spacing  is  such  that  mutual 

influence  becomes  important. 

(b)  The  rings  are  all  of  the  same  size  and  design. 


J 
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(c)  The  term  "a"  will  be  taken  to  be  the  outside  of  the 
shell  instead  of  mean  radius.  For  shells  with  a 
large  radius  to  thickness  ratio,  the  effect  is 
negligible. 

Boundary  Conditions:  (a)  The  radial  deflections  of  the  shell  and  ring  are 

equal  at  the  junction. 

(b)  The  slope  of  the  shell  adjacent  to  the  ring  is  zero. 

Analysis:  As  a  consequence  of  the  proximity  of  the  ring  stiffeners,  the  long 
cylinder  must  be  treated  as  composed  of  a  number  of  short 
cylinders,  each  of  which  has  a  length  i.  The  discontinuity  forces 
and  moments  at  the  junction  of  the  short  circular  cylindrical  shell 
and  the  circumferential  ring  stiffener  may  be  taken  to  be  identical 
to  those  shown  in  Fig.  B.3.2.  In  this  case,  however,  the  deflec¬ 
tion  of  the  short  circular  cylindrical  shell  at  the  junction  produced 

by  the  forces  p,  Q  and  M  is,  from  Reference  4, 
o  o 


The  radial  displacement  of  the  circumferential  ring  stiffener  has  been 
calculated  previously  and  is 


In  accordance  with  boundary  condition  (a),  we  must  have 


.  y  »•■>«  v  j»*, 

E^>  \s.’»upC  +  sU|?*  /  EW»  \s^l,pi+supl  /  im'  ' 

»  (aa.+pcWbVa*  ] 

cE  Vbx-ax  / 

(B.4.1) 


The  slope  of  the  short  circular  cylindrical  shell  adjacent  to  the  ring  produced 
by  the  forces  p,  Qq  and  may  also  be  obtained  from  Reference  4.  Specifi¬ 
cally,  it  is 


2  Softly  stnkfU-sUfjt  \  4M(  /  CotUfl-  \ 

EU  V  4  \  S»KU^f+Smj|4  / 


This,  according  to  boundary  condition  (b)  should  be  zero.  Hence 


a.- 


(B.4.2) 


If  we  now  substitute  Eq.  (B.4.2)  into  Eq.  (B.4.  1),  an  explicit  solution  for  the 
circumferential  bending  moment  Mq  will  be  obtained.  Indeed, 


[«■ 

('-!)- v>l 

V 

J  .  (B.4.31 

k&h 

A-46 


whore 


Fig,  B.5.1  -  Long  Circular  Cylindrical  Shell  with  a  Flat  Head  Closure 
under  Internal  Pressure.  The  forces  and  moments  are 
shown  in  the  positive  sense  above. 

The  following  compatibility  equations  between  cylinder  and  head  are  obtained 
from  Reference  5,  and  are  modified  to  conform  to  the  notations  of  Fig.  B.  5.  1. 


—  a,+  ar  +  = 

*  5  6  4 


i  q  r-  |  V 

- -a.+ - a*  (B-5-l) 
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Short  Circular  Cylindrical  Shell  with  Equal  Thickness  Flat 
Head  Closures 


Fig.  B.6.  1  -  Short  Circular  Cylindrical  Shell  with  Equal  Thickness 
Flat  Head  Closures  under  Internal  Pressure 

The  signs  convention  and  notations  of  Fig.  B.6.  1  will  be  used  in  the  present 
analysis.  For  the  case  of  a  short  cylindrical  shell,  the  radial  displacement 
(positive  inward)  6,  and  the  rotation  -6'  at  the  junction  due  to  the  forces  p, 
Qq  and  Mq  are,  respectively, 
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The  dimensionless  form  of  these  two  equations  may  be  chosen  to  be 


a4  + 


M,  - 


4f>«.* 


Qs  + 


a°  a. 


7  r><\ 


v  _  _M;.r  a-  T 

8|9<»3pZ^,  4j9At  4  5 


(B.6. 1) 


(B.6.2) 


where 


r-  __  _  J  j'j^Y  \ 

b  Z  \  ^»|/  \  sisU  ^  l  4  / 


(B.6.3) 


These  "influence  numbers”  may  be  substituted  in  lieu  of  the  influence 
numbers  a^,  a^,  a^,  b^  and  b^  for  the  long  cylindrical  shell  in  Eqs.  (B.5.4) 
and  (B.  5.  5)  to  give  the  expressions  for  the  shear  force  and  bending  moment  at 
the  junction  where  the  flat  heads  abutt  the  cylindrical  shell,  respectively. 


Specifically,  they  arc 


M,=  4f>«1 


(q4~  °s)Qy  W6)  +  Qj) 

_  •“  ®  tX k4” ^5) “  C° i“aeX'Bt“  W). 


C^3~  b«.)  *  W(g5~q  1) 

(<V5SX  bz-w5)-  (<V«6X  b  -  b4) 

Defining 

CotU  (ll  +  C 
s  *  *  W  4  S'*' 
cosUpC  -c+*pZ 
S.*v,(j<£+sup* 


and  substituting  Eqs.  (B. 4. 4)  and  (B.5.3)  into  Eqs.  (B.  6.3),  we 
following  relationships: 


“4  = 

a4 

a*K 

*G  = 

f' 

II 

<T| 

W* 

II 

(B. 6.4) 


(B.6.S) 


(B.  6.6) 


obtain  the 


(B.6.7) 
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Equations  (B.6.7)  express  the  correspondence  of  Lhe  influence  numbers 
between  the  long  and  short  circular  cylindrical  shells.  Substitution  of  these 
relationships  into  Eqs.  (B, 6. 4)  and  (B.  6.  5)  yields 


~  (<V  M»)  ♦  _ 


(B.6.8) 


(as~  +  ^3 ^t) 

(a  f  aM  K-  -  (fit-  aADCK-vW . 


7.  Short  Circular  Cylindrical  Shell  with  Unequal  Thickness 
Flat  Head  Closures 


>C*o  :*»/ 


Fig.  B.7. 1  -  Short  Circular  Cylindrical  Shell  with  Unequal  Thickness 
Flat  Head  Closures  Under  Internal  Pressure 


For  the  short  cylindrical  shell,  the  radial  displacement  (positive  inward)  6, 
and  the  rotation  -6*  at  the  junction  where  x  0  due  to  the*  forces  p,  Qj,  0^, 

M y  and  are,  respectively, 

Gl  v  (sUV^JE  -  S»*rv3ftc»J  pt)  +  Q2(  sUUpl  CQlU  2t) 

,‘mo  2(a*o(sUk*^fi  -  s:»'tff) 

M | ( si «ejB C  +  s« »l p sjrtU jSC 

2  |3Z  O  (s>  »V,zp{  _  si  ) 

s'  —  2Qzs;*pls;A(U  -  Qt  4  s:^w*pg) 

***  2  pzo  (s;.W*p  &  -  s  :«zpC) 

Mt  (s.»fU  c»apl  +  sn.v>p<  c»»UfiC)  -  Mt(si>->pCo»sV»p< •»  c«,ft  »;„upO 
pD(si«Wlpi- (B.7.2) 


ZEk 


(2**0 


(B.7.1) 


For  the  junction  x  =  i ,  the  radial  displacement  and  the  rotation  are, 
respectively, 


si4^t- sup#  udkj(t)4  az(&^Uplc»^{t<-sUf4o«->|t) 

S*"*  2^So(sJ~W*fi-  »:**pA) 

+  7^tftt  f  2Mi  liifi  *"ufl  ~  M*(s:~wfo  ♦  *■"¥*)  ]  _  r?l(z_y) 

E^[  x^u'ft  -  s.v‘pl  J  2EU,v 


(B.7.3) 


4p*Q*T  m,  »  c»^ts.vwpl)-»A2(s»^ti>«>(rt'»».vv^<c«»Vfl) 

(B.7.4) 
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Equations  (B.7.1),  (B.7.2),  (B.7.3),  and  (B.  7. 4)  may  bo  arranged  in  the 
following  dimensionless  forms: 


J 


Let  us  now  define  the  following  parameters: 

_  sinh^  pi  +  sin^  pi 

°20  =  .  2  .1 

sinh  pi  -  sin  pi 

„  sin  pi  sinh  Pi 

4  sinh4  pi  -  sin  pi 
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p  _  ginh  p I  coBh  PI  -  gin  PI  con  6 I 
sinh2  pi  -  sin2  pi 

o  =  £2f  P j  pi'  -  gin  g|  cosh  B1 

18  2 - T - c— 

sinh  pf  -  sin  pi 

p  =  sin  P*  cog  pi  +  sinh  pi  cosh  p I 
sinh2  pi  -  sin2  pi 

p  =  sin  pi  cosh  pi  +  cos  PI  sinh  01 

19  sinh2  pi  -  sin2  P I 

Furthermore,  let 


(B*  7.  9) 


(cont. ) 


j 


j 
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(B.7.10) 


• 


Aio' 


Equations  (B.  7.  9)  and  (B.  7.  10)  may  now  be  substituted  into  Eqs.  (B.  7.  5) 
through  (B.7.8)  to  give  the  following  set  of  dimensionless  equations: 


4px  **  4  ]»a 

ip: 


(B.7.11) 


A,t) 


(B. 7. 12) 


6 


(B.7.13) 


8pVWj» 


4^  +  4^*»  (  i?)  +  ^  (" 2AkA*>) 


+  2«^  K,A*-'1 


(B.7.14) 


From  Reference  5,  we  find  the  following  dimensionless  equations  for  the 
heads: 


Q, 

Jii.p 

4px  “ 

Zj>A 

_ 

Mz. 

— =:  4 

4^»at 

Z-J9A 

_ 

a- 

8  1» 

4f>*  ' 

EU2  *' 

♦  -Si- 

8  |EV'Uj|i> 

V* 

2jxv 

T  •“‘zi 


(B.7. 15) 


(B.7. 16) 


(B.7. 17) 


(B.7. 18) 


where 


=  2(l-y) 


(1-2/) 


*»7  ~ 


gQ-*0 

plV>,Uz 
.  50->) 
zpVhj 


=  fc, 


(cont.) 
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^iZ 


\6|S 

&0-y) 

.  30-y) 


(B.7.19) 


By  equating  Eqa.  <B.7.11)with  (B.7. 15),  (B.  7.  12)  with  (B.  7.  16),  (B.7.  13) 
with  (B.7,  17)  and  (B.7.  14)  with  (B.7.  18),  we  arrive  at  the  following  set  of 
equations  in  Mj »  M^,  Qj#  andQ^.  Specifically,  they  are: 

AOA:  nO - 


^j^(KAfS)+  *» 


(B.  7.20) 


Equation  (B.  7.  20)  may  also  be  conveniently  expressed  in  matrix  form  by 
defining 

(Mil  (Ab-AjI 

l„l  M2  .  U  JA16-A4I 


%  ;  |v|  -  P‘*J 

Q,  U  I  A,q 


(conU) 


o 
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Fig.  B.8*  1  -  Forces  and  Moments  at  the  Junction  of  the  Long  Circular 
Cylindrical  Shell  and  the. Ellipsoidal  Head  Closure  under 
Internal  Pressure 


A-60 


Fig.  B.  8. 2  -  Signs  Convention  for  Middle  Surface  Radial  Deflection 
and  Rotation  at  the  Junction 

The  radial  deflection  and  the  rotation  of  the  long  circular  cylindrical  shell  at 
the  junction  produced  by  the  loads  p,  Qq,  and  are,  respectively 


(B.8.1) 
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A- 62 


from  which 


(ft>'  ?*0v^+p.>«)a*  ■  t*] 


(B.8.8) 


^ftP.+  P.OM‘+  2^>«'P.*D.)Q*  °  (B.8.?) 


Simultaneous  solution  of  Eqs.  (B.  8. 8)  and  (B.  8.9)  yields 


(B.8.10) 


(B.8.11) 
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These  formulas  for  Mq  and  may  be  put  into  the  following  dimensionless 
forms: 


Short  Circular  Cylindrical  Shell  with  Equal  Thickness 


Ellipsoidal  Head  Closures 
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The  signs  convention  and  notations  of  Fig.  B,  8.  2  will  be  used  in  the  present 
analysis.  For  the  short  circular  cylindrical  shell,  the  dimensionless  forms 
of  the  equations  for  radial  displacement  6^  and  the  rotation  6^  at  the  junction 
due  to  the  forces  p,  Q^,  and  are  given  by  Eqs.  (B.  6. 1)  and  (B.6.2), 
respectively.  These  are  listed  in  the  following  forms: 


By** 

tftt1 


<V 


<a„ 

2  ha 


«4 


3 


(B.9.1) 


(B.9.2) 


where  a^,  a^.  a^,  and  are  a*  defined  in  Eq.  (B.  5.  3);  by  Eq.  (B.4.4) 
and  and  by  Eq.  (B.  6. 6). 

For  the  ellipsoidal  head  closure,  the  radial  displacement  6  and  the  rotation 

e 

'w  6*  at  the  junction  due  to  the  forces  p,  Q  ,  and  M  are  given  by  Eqs.  (B.8.4) 

e  o  o 

and  (B.  8.  5).  These  two  equations  may  be  written  in  the  following  dimension¬ 
less  forms: 


tu,se  M0  <ao  ^ 

- —  - -  a  ♦ — —  az  +  a3 

4-f>«.  2jj« 

-  ihl  s«  M.  .  Q.  , 

4  pa*  1  zt*  z 


where 


(B.9.3) 


(B.  9.4) 


(cont.) 
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(B.9.5) 


In  accordance  with  the  strain  compatibility  equations  (B.8.7),  we  have  the 
following  equalities: 


J 


•—  ar>  +  a  £i  * 
4  4jx*  J  15  zpn  6  4 

Aya1  4  5  2j»a  » 


M. 


£2,+ 


4f«*  * 

-^b,+ 

4pa  1 


a. 

2jya 

a, 

2p6. 


Oj_4  rt; 


J 


from  which 


M„=  4p©1 


( a4-  a3)(  \o7- bs&l3) 


(W^JOv^b)-  V  b4SLl3) 


(B.9.6) 


a0=  ^ 


(a  a  -  b,  - 

(a a  5&iX  bj~  b^j)  -  (p2-  b,“ 


(B.9.7) 


J 
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10.  Short  Circular  Cylindrical  Shell  with  Unequal  Thickness 
Ellipsoidal  Head  Closures 


a,  a,  a2 


Fig.  B.  10.  1  -  Short  Circular  Cylindrical  Shell  with  Unequal  Thickness 
Ellipsoidal  Head  Closures  under  Internal  Pressure 

The  signs  convention  of  Fig.  B.  8.2  will  be  adhered  to  in  this  analysis.  For 
the  short  circular  cylindrical  shell,  the  radial  displacements  and  rotations 
at  the  junctions  where  x  =  0  and  x  =  I,  respectively,  due  to  the  forces  p,  Q^, 
Q^,  Mj,  and  are  given  by  Eqs.  (B.7.  11)  through  (B.7.  14)  in  dimensionless 
forms  and  will  not  be  repeated  here. 

For  the  ellipsoidal  heads,  the  corresponding  dimensionless  forms  of  the 
radial  displacements  and  rotations  are: 


_  = 


2-pA 


(B.  10. 1) 


•c\.  S 
*.=  £ 


2r 


(B.  10.2) 
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J 

(B.  10.  3) 


(B.  10.4) 


J 


(cont.) 

J 


o.  =- 


= 


ibl 


,2.0-y*) 


zw. 


\2  0-^) 


(B. 10.5) 


Equating  now  Eqs.  (B. 7.  1 1)  with  (B.  10.  1),  (B.7.  12)  with  (B.  10.2),  (B.  7.13) 
with  (B.  10.3),  and  (B.7. 14)  with  (B.  10.4)  in  accordance  with  the  strain  com¬ 
patibility  equations  (B.  8.7),  we  arrive  at  the  following  set  of  equations  in  Mp 
M^,  Op  and  Q^.  Specifically,  they  are: 

+  +  *•  *o-*» 


£r( riCMtf ££r0 £tl**m*J  -  ^-^4 


4r 


zr 


zr 


^rCNV  -  ° 


Let  us  now  define  the  following  matrices: 


M 


a. 


W- 


*k'N- 

o 

O 


(cont.) 
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J 


4j>a1 

4|»a4 

4y*»l 

4| >«* 

4rl 

4^>al 


^3$**^*  ^3-0* »» 


2|>a 

2|>A 

^^•6" 

2pa 

2|»A 

zA^biXrt 

Z 

2r 

Zj^  ZjjA. 


(B.  10.6) 


The  coefficients  Aj  through  Ajq  in  the  above  equations  are  defined  in 
Eq.  (B.7.  10)  where  it  is  understood  that  p  =  P^.  The  functions  through 
n^g  are  identically  defined  in  Eq.  (B.7.9). 


Circular  Cylindrical  Shell  with  a  Hemispherical 


Head  Closure 


a,  a0 


Fig.  B.  1  L  1  -  Forces  and  Moments  at  the  Junction  of  the  Long  Circular 
Cylindrical  Shell  and  the  Hemispherical  Head  Closure 
under  Internal  Pressure 


Fig.  B.  11.2  -  Signs  Convention  for  Middle  Surface  Radial  Deflection 
and  Rotation  at  the  Junction 
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The  discontinuity  forces  and  moments  at  the  junction  of  the  long 
circular  cylindrical  shell  and  the  hemispherical  head  closure  under  the 
action  of  internal  pressure  can  be  readily  obtained  from  Eqs.  (B.8.  1Z)  and 
(B.  8.  13)  developed  for  the  case  with  the  ellipsoidal  head  configuration;  it 
being  necessary  to  replace  the  quantity  a^/b^  by  unity.  In  so  doing,  we  find 


(B.  11.2) 


1Z.  Short  Circular  Cylindrical  Shell  with  Equal  Thickness 


Hemispherical  Head  Closures 


Fig.  B.  1Z.  I  -  Short  Circular  Cylindrical  Shell  with  Equal  Thickness 
Hemispherical  Head  Closures  under  Internal  Pressure 


o 


The  formulas  developed  for  the  casi*  of  the  short  circular  cylindrical  shell 
with  equal  thickness  ellipsoidal  head  closures  arc  readily  adaptable  to  the 
present  configuration  by  substituting  unity  in  place  of  the  quantity  a^/b^. 

For  example,  the  quantity  in  Eq.  (B,  9.6)  for  the  bending  moment  MQ  and 
in  Eq.  (B.9.7)  for  the  shear  force  must  now  be  replaced  by 

a3  =  -|(l-v)  (B.12.1) 

13.  Short  Circular  Cylindrical  Shell  with  Unequal  Thickness 
Hemispherical  Head  Closures 


o 


fl,  a,  a,  a. 


Fig.  B.  13.  L  -  Short  Circular  Cylindrical  Shell  with  Unequal  Thickness 
Hemispherical  Head  Closures  under  Internal  Pressure 

Formulas  (B.  10.6)  and  (B.  10.7)  for  the  determination  of  the  bending  moments 
Mj  and  and  the  shear  forces  and  developed  for  the  case  of  a  short 
circular  cylindrical  shell  with  unequal  thickness  ellipsoidal  head  closures  are 
also  applicable  to  the  present  configuration  with  the  exception  that  the  quanti¬ 
ties  Ajj  and  defined  in  Eq.  (B.  10.5)  must  now  be  replaced  by 

A13  =  Al6=  -|{  1  -  V)  (B.13.1) 
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whore 


For  the  conical  head,  the  corresponding  dimensionless  forms  arer 


A'jxx1 

4^a«-  a' 1 

2^ 

"^Tb.+ 

Q. 

CD 

I- 

X 

' 

2-jo<v 

V»s 


(B.  14.7) 


(B.  14. 8) 


The  quantities  a^,  a^,  a^,  bj,  b^,  and  b^  are  dimensionless  "influence 
numbers"  for  the  conical  head  and  can  be  obtained  from  either  References  Z 
or  3.  In  the  notations  of  Fig.  B.  14.  1  and  the  signs  convention  of  Fig.  B.  14«2, 
these  influence  numbers  become 


2,  ~~ 


(cont.) 
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S»«  o< 


az  = 


«»- 


**& 

-  ^T\  fl 

C+2*><* 

nr* s<crt  *  s<c«x 


b  -  - 


26  f.  Wot 

C+2vG  4 


_  k/2, 

C+2*>4 


-  CQ  +  »)i,  ^  hiiinsl  .3  S*c<*e«s«c» 

c  4  2  ?: 


(B.  14.  9) 


where 


S.(b«rt'  S.  b«.'t?„-  be/ 

0~*  t.)\  (bet'  ?.)* 

$ .(b*,r*  $J®*fI**  *  *•) 

(b^f+O^O* 

^ ■^C«tol  C  •$«&<* 


(B.  14. 10) 


For  strain  compatibility,  we  must  have 


6  =  -6 
c  co 


6'  =  6’ 

C  CO 


(B.  14.11) 
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Hence, 


M.  Q, 

V  °*+  <v‘,‘ 

4f»*  4  2^  * 


M.  „  6». 

- -fl, - g  -  A, 

2fA  1 

M,  ,  Q 

4^  '  *pA  *4  5  <B.  14. 12) 


(a,i.a4)t 

0. 

4jsa*- 

"2  jo<X 

<3U 

4-p^ 

2]»a 

OV«s)  «  -(«»♦«,) 


from  which 


(B.  14. 13) 


M.m  4j.*1  (a,<°c)(^-bs)- b.(a,va5) 

L(«,«*xts-io+ (k.rt.4xa«*«s)J  (B-U' 

a.-  *r«  _S(ai-^>«)->(att«t)(b4.-b.)  1  m 

L  (<*i*a«Xbs'*>0  +  (b,-b1)(att  «s) 

Th.  Bee.el-Kelvin  function,  bur,  bei,  ker,  kei.  and  th.ir  f;„t  derivative, 
hau.  been  ..ten.iv.ly  tabulated  by  Lowell  (Ref.  7).  Higher  order  B....1- 

Kelvin  function.  ...  related  lh.  aforementioned  function,  by  th.  following 
recurrent  formulas  (Ref.  8): 

b*’r'1  “  V*'  b*i’*) 


(cont.) 
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ual  Thickness  Head 


Closures  of  Ellipsoidal  and  Conical  Shape 


Q, 


kb- 


T 

a 


M,  M, 


XSO 


M,  M, 


»  « 

a 

1  F* 

V’ 

1 1,1 1 

FxF 
1 1 1 

r  V 

> 

a.  a* 


Fig.  B.  15.  1  -  Short  Circular  Cylindrical  Shell  with  Unequal  Thickness 
Head  Closures  of  Ellipsoidal  and  Conical  Shape  under 
Internal  Pressure 

The  signs  convention  of  Figs.  B.8.2  and  B  14.2  will  be  adhered  to  in  this 
analysis.  For  the  ellipsoidal  head,  the  dimensionless  forms  of  the  radial 
displacement  and  rotation  have  been  derived  earlier,  and  are: 


4 1*  a1 


4|>a* 


kit* 


M. 

4|»a* 


^«7  * 


r 


(B.  15.  1) 


(B.  15.2) 
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whore 


"  v,l  V3 C'_jyl) 

V  "  ITV/30-’'*) 

n«  i 

Ns~  -sC2--^-^) 
'“'n =  ^7  V’sO-f*)' 


N»- 

v  ^  ; 

zpelaUt 

f«  - 

_ 

4«ln>c 

a 

1*0 

4a1o, 

2  ,  Z 

a 

EVt 

o>  = 

c 

iz(\-ya) 

tv,? 

— 

For  the  short  circular  cylindrical  shell,  the  dimensionless  forms  of  the  radial 
displacements  and  rotations  at  x  =  0  and  x  =  f  due  to  the  forces  p,  Qj ,  Q^, 

Mj,  and  are,  from  Eqs.  (B. 6.  11 )  through  (B.  6.  14): 


Q» 


(B.  15.4) 


(cont.) 


I 


Finally,  for  the  conical  head,  the  dimensionless  forms  of  the  radial  displace¬ 
ment  and  rotation  at  the  base/of  the  cone  due  to  the  forces  p,  Q^»  and  are 
represented  by: 


u 


=  _  Mi_a  - 
4f>a.1  4  pa1  1  zpa 


(B. 15.9) 
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t 

'  and 


M 


HZI 

M.» 

4j»^t 

'  2.  j»ft* 

2  j*a 

2^0. 

^7-^(3 

4-p«l 

4^ 

2  pa 

fa 

A, 

4pftT’ 

2pa 

2^a 

7 1, 

k|%&i»+W 

4p*»l 

fa 

2. 

(B.  15. 12] 


The  four  simultaneous  equations  listed  above  may  be  written  in  the  form 


V 


from  which 


(B.  15.13) 


■nJ 
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16.  Short  Circular  Cylindrical  Shell  with  Unequal  Thickness  Head 
Closures  of  Hemispherical  and  Conical  Shape 

If  the  ellipsoidal  head  of  Fig.  B.  15. 1  is  replaced  by  a  hemi- 
spherica  head  of  radius  a  and  thickness  hj,  the  discontinuity  forces  Q^, 

Mj  and  associated  with  this  new  configuration  can  be  readily  calculated 
by  the  formulas  developed  for  Case  15.  In  this  case,  it  is  merely  necessary 
to  replace  the  quantity  defined  in  Eq.  (B.  15.  3)  by 

(B.16.1) 
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V 
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